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Abstract 

We study the tunnelling trough a potential barrier of the system of two 
quantum correlated particles. The system is considered in one dimension. 
The interaction with the barrier and between particles is approximated by 6- 
potentials. Assuming that the particles have the same masses, we reduce the 
problem to the set of independent planar scattering problems corresponding 
to two-body flux symmetry with respect to the barrier. In order to solve the 
problems we apply the method of Sommerfeld-Maluzhinets integral transfor- 
mation, which requires the solution of a set of functional equations. We show 
that the sub-problem which is antisymmetric with respect to the center of 
mass allows the solution in frames of the Bethe ansatz. So we give the exact 
solution of this problem. For the subproblem which is symmetric we describe 
the class of solutions for the functional equations and give an approach for 
obtaining the solution. 

1 Introduction. 

The properties of low-dimensional semiconductor structures have attracted much 
attention in the past few years both experimentally and theoretically. Being the 
quantum objects by nature the ID and 2D nanostructures can be described by an 
appropriate Schrodinger equation with the suitable boundary conditions. There exit 
a huge literature on the subject however the most results obtained deal with the 
solution of 1-body or 2-body Schrodinger equation and with the physical quantities 
that can be described in terms of these solutions. To our knowledge there are only 
few results concerning the 3-body systems in nanostructures while there are a lot of 
physical phenomena which need the solution of 3-body quantum problem. Let us 
mention only two of these phenomena: 1) Observation of anomalies in the optical 
spectra of quantum wells have been attributed to the formation of negatively-charged 
excitons [23] . Such complexes can arise when exciton E captures one electron to 



form a negatively-charged exciton E + e — > E~ |21| . In ref. [21] the properties of 
such 3-body complexes in ID quantum wires have been studied by means of exact 
solutions of ID 3-body Schrodinger equation with the Calogero type interaction 
|22j . 2) Since the discovery of the vanishing transmissivity of a tunnel barrier in 
a ID quantum wire due to the electron-electron interaction [22] > new interest has 
emerged in the transport properties of ID electron systems. Indeed, the influence of 
the electron correlations shows up strikingly in non-linear current-voltage relations 
which are investigated experimentally in narrow semiconducting wires [21]. 

From the mathematical point of view the type 2) of physical phenomena can be 
described by 3-body Schrodinger equation in one dimension if one treats the local 
electron interaction via 5-like potential and the quantum barrier (well) as the third 
infinitely heavy particle which interacts with the electrons by ^-function potential 
as well. Solving the appropriate 3-body ID scattering problem one can calculate in 
terms of scattering amplitudes such properties of tunneling barrier in quantum wire 
as conductance, current-voltage characteristics, etc. This scattering problem can be 
solved either numerically or analytically i.e. exact. The present paper is devoted 
to the exact solution of the problem mentioned above by means of the reduction of 
the original 3-body Schrodinger equation to some Riemann-Hilbert problem. Using 
the symmetries of the system we split the latter into two subproblems which are 
either symmetric with respect to the center of mass or antisymmetric. For the 
antisymmetric case we construct an exact solution in terms of the Bethe anzatz. 
For the sub-problem which is symmetric we suggest a new approach which allows 
to describe the functional classes of the solutions and reduce the problem to some 
generalized Dubrovin's equations. 

The paper is organized as follows. In Sec. 2 we briefly describe the mathemati- 
cal model for the correlated tunneling of two electrons through a barrier. In Sec. 3 
by means of Sommerfeld-Maluzhinetz transformation we construct the matrix func- 
tional equations which have to be solved. In Sec. 4 we give an exact solution of 
antisymmetric sub-problem in terms of the Bethe anzatz. Sec. 5 is devoted to the 
study of the Riemann-Hilbert problem for the symmetric subproblem. Sec. 6 contains 
some concluding remarks. 

2 The Model 

The model we study here has been suggested by K. Lipszyc [7j in 1980. However, 
regardless the simple way the problem is posed, no closed-form solution is known 
until now. Here we shall come closer to the exact solution of the complete prob- 
lem and present the solution for one of two sub-problems it is reduced to. Since 
the problem was carefully described in [7j, we shall allow us to abandon the de- 
tails and present here just the topics essential for understanding of the following 
considerations (terms, notations, etc.). 

We consider the system of two quantum particles with the same mass in one 
dimension interacting via delta-function potential with the strength 7. In addition, 
each particle is interacting in the same way with the fixed delta-function barrier 



with the strength T. Let us denote the coordinates of the particles with respect 
to the barrier by r% and r 2 , their masses by mi and m 2 , choose the units so that 
mi = m 2 = h/2 = 1. The Srodinger equation for the wave function u of the system 
can be written in the form: 



-A - T<5(ri) - T5(r 2 ) - 7<5(ri - r 2 )]u(r 1 , r 2 ) = Xu(r 1 , r 2 ). 



(1) 



The correct mathematical formulation of this equation and the corresponding scat- 
tering problem is given in terms of extension theory fSJ Ej and was described for 
similar problems in ^HlEl El- Here we shall come directly to the equivalent bound- 
ary value problem. To do this, let us introduce the polar coordinate system for 
(ri,r 2 ) = (r cos(6 l ), r sin(6 1 )), denote the function u in polar coordinates by u and 
define 

Lf](r,0) :=/(r,0 + O)-/(r,0-O). 
In this terms, the boundary value problem we shall study looks as pQ : 



i-u 



—Au = Xu 
= — rTu(r, <fi) 



</>=0,7r/2,7r,37r/2 

2r7«(r, 6) 

y ^' 0=7r/4,57r/4 



(2) 



Taking into account the symmetry of the system we split the solution of (0) in 
two parts, symmetric and anti-symmetric with respect to the particle interchange: 



w(ri,r 2 ) = M_(ri,r 2 ) +u + (r 1 ,r 2 ) 
u-(n,r 2 ) = -w_(r 2) ri) = (u(r 1 ,r 2 ) - «(r 2 ,ri))/2 
u+(ri, r 2 ) = +u + (r 2 , n) = {u(ri, r 2 ) + u(r 2 , n))/2 



(3) 



The solution of Q for is trivial since the role of the second boundary condition 
(the interaction between particles) vanishes being replaced by the Dirichlet boundary 
condition u(x, x) = u(r, vr/4) = 0. The expression for U- can be found as: 



w_(n,r 2 ) = /i(ri)/ 2 (r 2 ) - A(r 2 )/ 2 (ri) 



(4) 



where j\, 2 are the solutions of the one particle problem with 5-interaction . 

We shall concentrate on the solution of the boundary problem (0) for u + . Due to 
the symmetry of u + the boundary value problem (j2J) can be expressed simpler in the 
centre of mass frame: x — (r\ + r 2 )/y / 2; y — (r\ — r 2 )/^/2, where x is the v^2 times 
coordinate of the centre of mass and y is the distance between particles divided by 
a/2. Unless otherwise stated we shall use this coordinate system for u + , the function 
u + will represent the same function in polar coordinates. In this notations, one has: 



9$ 



-Au + = \u + 
= —rTu + (r, q 

) = —r^u + (r, 

0e (o,tt) 



$=7r/4,37r/4 



0=O,7T 



(5) 



The problem can be further simplified using its symmetry with respect to 
the transformation x — > — x (0 — > tt — 0). Let us split the function w + accordingly: 



u+{x,y) = v-(x,y) + v+(x,y) 
v-{x,y) = -v-(-x,y) = (u+(x,y) - u+(-x,y))/2 
v+(x, y) = +v + (-x, y) = (u+(x, y) + u+(-x, y))/2. 



(6) 



The representation © splits (0) into two independent problems for v±: 



-^v±\ (r,7r/4) = -rTv±(r,n/4) 
■§-v± + (r, 0) = —r r yv±(r, 0) 
l)i;_(r,7r/2) = 0; 2) ^ + (r,vr/2)=0 

0e (o,tt/2) 



(7) 



These are the boundary value problems we shall deal with in the rest of this paper. 
In order to be rigorous we must add that we shall investigate the corresponding 
scattering problems. That is we shall look for the solutions of (J7J) being the combi- 
nation of one incoming plane wave e~ lkrcoB ^~M , O £ (0, vr), a set of outgoing plane 
waves, 0o ^ (0, 7r), and a function satisfying the standard scattering conditions at 
r — ► oo. The detailed description of how the scattering problem should be posed 
can be found in extensive literature on the subject [TJ IH 13 13 UH UHj- 

As a conclusion to this section, let us briefly discuss the physical meaning of the 
symmetries we applied. The best visualization can be achieved if we discuss the 
meaning of the corresponding incoming plane waves. The plane wave in the system 
(J2J) corresponds to the incoming, that is in which the particles are directed towards 
each other and towards the barrier, flow of pairs (since we consider the interaction 
in pairs, we should consider the pair of particles as one) of particles in which the 
momenta of the corresponding particles are the same. That is all particles of type 
1 have the momentum k\ and all particles of type 2 have the momentum fc 2 . We 
should mention that though the particles are assumed to have the same masses and 
interact with the barrier in the same way, they can, in principle, differ by some other 
properties (say by spin). The plane wave in the system (jSJ) corresponds to two fluxes 
which have the momenta or particles interchanged and have the amplitudes which 
are either the same or differ by —1. In the case of absolutely identical particles, this 
would correspond to their statistics. In the last case, the plane wave in the system 
(0) corresponds to four fluxes coming from two different sides to the barrier, by two 
from each side, and whose amplitudes are either the same or differ by —1. From the 
physical point of view, this situation as separate seems to be hardly experimentally 
realizable but may also give a better understanding of the process. In this paper 
we give the exact solution of the scattering problem (J7|) with boundary condition 1) 
and approach the solution of the (JJJ) with condition 2). 



3 Sommerfeld-Maluzhinetz Integral Transforma- 
tion and The Functional Equations. 



To discuss the solutions for (JJJ), K.Lipszyc ,7] used the method of Sommerfeld- 
Maluzhinetz (SM) Integral Transformation ^U|. The Maluzhinets type functional 
equations were derived, however the solution was said to require some new unknown 
functional class to be used. Thus the problem was closed until further progress in 
functional analysis. Here we shall also use the SM transformation, however we shall 
study the functional equations in the form of McGuire [21 E] . This way of using the 
SM transformation was effectively applied to several scattering problems El ESj • 
In our previous papers [T31 C] we proved several general statements which described 
the functional classes needed to solve some SM functional equations (the McGuire 
type) and the way one can construct their solutions. Here we apply these statements 
to the concrete problem in question. 

Let us apply the SM transformation to the function v(r, 0) considered on two 
domains: v(r, 0) = v(r, (ft), G (0, 7r/4) and v(r, 0) = v(r, + vr/4), 6 (0, 7r/4). 



v (r, 0) = J c e~ ikr cos a (g x (a + 0) + g 2 (a + vr/4 - 0) ) da 
v(r, 0) = f c e -* krcosa (g 3 (a + 0) + g A (a + tt/4 - 0)) da 



where c stands for the standard Sommerfeld contour (TU1 El EE] • Using the represen- 
tation (JHJ) and the Maluzhinets theorem JU] the boundary value problem (J7J) can be 
transformed into a set of functional equations: 

iksm(a) (g 3 (a) - g±(a + 7r/4) - g x [a + 7r/4) + g 2 (a)) = -T (g 3 (a) + g 4 (a + tt/4)) 
93(a) + g 4 (a + 7r/4) = g x (a + tt/4) + g 2 (a) 
iksi\i{a) [gi{a) - g 2 (a + vr/4)) = -7 (g x (a) + g 2 (a + vr/4)) 
5f 3± (a + vr/4) ^ g 2± (a) = 

(9) 

Here we add the + or — signs to the functions corresponding to v + and v _ respec- 
tively. 

If we introduce the notations g(a) = (g x (a), (74(a), gz{a), g^a)) 1 the functional 
equations can be rewritten as 

g(a + tt/4) = M(a)g(a) (10) 

where the meromorphic matrix function M is given by 



M±(a) 



/ n n 2ik sin(q) T 

' 2ifcsin(a)-r 2ifcsin(a)-r 

pi pi T 2ik sin(a) 

U U 2ifcsin(a)-r 2ifcsin(a)-r 

±1 



ifc sin(ct) 



v 00 1 



ID 



To study the functional equations (fTUj) we need some algebraic results obtained 
in our previous papers PIE] for the generalized class of similar functional equations. 
Below we apply these results to the problem under consideration. Namely, let M (x) 



and F(x) be meromorphic matrices depending on the complex variable x and such 
that for some n G N 



M (x + n) = M (x); F(x + n) = F(x). (12) 

The functional equation (fTUj) with the argument shift 7r/4 can be reduced to the equa- 
tion with the argument shift 1 and the matrix coefficient having the period n instead 
of 27r. The properties of matrix coefficient M(a) from (jllj) can be reconstructed if 
one poses n = 8, x = 4a /n. We have considered the matrix finite-difference equation 

g(x + 1) = M{x)g{x) + F (x) . (13) 

We introduced the notation for ordered product of some matrices {Mj}, namely 
T Ui=i Mi = M k M k _ 1 ...M l for k > I and | Ui=i M { = I for k < I. In this terms the 
n-th iteration of the equation (JT5j) reads as 

g(x + n) = N(x)g (x) + E (x) (14) 

where we denoted N (x) =| 117=0 M {x + i) and S (x) = E7=o T U]=!+i M (x + j) F(x+ 
i). Due to the definition of iV one has 

N (x + n) = N (x) ; E (x + n) = E (x) (15) 

N(x + 1) = M (x) N (x) M~ l (x) . (16) 

So the eigenvalues of N (x) are of the period 1. Then we considered two special 
cases! a ) iV (x) = I b) det (N (x) — I) ^ which covered the three body problem 
with zero-range interactions and are suitable for the problem we consider. We need 
the following statements [T| I13j: 

Lemma 1 The general solution of inhomogeneous functional equation P^l) cor- 
responding to the cases a) and b) looks like 



a ) 9 \x) = %E (x) + g (x 

b) g(x) = (I-N(x))- 1 ~(x)+go(x) 



(17) 



where go (x) is a solution of homogeneous functional equation 

go(x + n)=N(x)g (x) (18) 



Using this statement we were able to proof the 

Theorem 1 The general solution of inhomogeneous functional equation (|J5y for 
the cases a) and b) looks like 



9 (x) = E?=o T npVi M (x + j) (&) F(x + i) + 



Er= T YTjZl M(x+ j) g (x + i) 



(19) 



h) g(x) = (I-N(x)) 1 ~(x) + 

) +E7=o 1 ^mzlM(x + j)go(x + z), 1 1 



where go (x) is a solution of homogeneous functional equation . For the case 
det N (x) ^ the solution can be also expressed as 



9 (x) = Er=o T Tip 1 M (x + j) (I-N(x + z))" 1 F(x + i)+ 

+ E»:o 1 tnSM(x + j) go (x + i). 1 ] 

The statements above show that the first step in the study of the equations (jHJ) 
should be the construction of the 8-th iteration of the matrix M(a). In order to 
simplify the calculations we introduce the notations 

T fc := T/k 

Ik ■= 27 A (22) 



and represent M as M = ^pj ^ where 



P + (z) 



±(z 2 - lk z - I) 

z 2 + lk z-l 



p(z) = Z 2 - 'JkZ - 1 

Q(z 



z 2 - 1 r fc ^) 



v2 



T k z z 2 - 1 
q(z) = z 2 - T k z - 1. 



(23) 



This allows us, after one iteration, to deal with the 2x2 matrices instead of 4 x 4 
ones, and consider the equation g(a + 2tt) = jV(a)g(a) as a set of two independent 
equations. 



4 Anti-symmetric case and the Bethe ansatz 

Iteration of the equations for t>_, i.e. the case which is center-mass anti-symmetric 
with respect to the barrier, shows a surprising result: 

AL(a) = /. 

As it was mentioned in pQ the latter result means that the Bethe ansatz can be 
used for the exact solution of the problem. One can construct this solution through 
either SM transformation (see e.g. jU]) or directly solving the finite linear system 
for the plane waves coefficients. Here we give the result of these calculations. 

We shall construct the solution as the decomposition of plane waves. Let us fix 
the wave vector (fci,/^) for one of the plane waves. Using the reflection rules we 
conclude that the complete set of the plane waves can be written as: 

<P e (0,tt/4) v4x,y) = D=i,..,8^e tk < x 

(7r/4,7r/2) v4x,y) = Ei=i,.., 8 B l e^ [ ) 



where 



1=1... ,8 





h, 


k2) ) 




-h, 


-k 2 ) 






-k 2 ) 


< \ 


-h, 


k 2 ) 




k 2 , 


h) 




-k 2 , 


-h) 




k 2 , 


-hi) 


V 


-k 2 , 


ki) ) 



(25) 



and x := (x, y). 

Let 3ftfci < 0, $lki < and $lki < 9ft/c 2 , that is the wave vectors (ki,k 2 ) and 
(k 2 , k 2 ) correspond to the incoming waves which represent two situations: first, the 
particles are coming from the same side of the barrier; second the particles are 
coming from different sides of the barrier. To get the scattering amplitudes one 
should solve the boundary value problem (JJJ) for A\ = 1; B5 = (first case) and 
A± = 0; B$ = 1 (second case). The results we obtained for the scattering amplitudes 
are given by equations (|2Tj) and (j27|). 



At A 1 = l; B 5 = we have: 



A 2 
A 3 
A, 
A 5 
A e 

A 7 
A 8 
B x 
B 2 

B 3 

£4 
Be 

B 7 
B 8 



2 7 fci 2 (7-H fc 2 )+2fci 3 (-i 7 +fc 2 )+fci (-i 7 +fc 2 ) (r 2 -2 fc 2 2 ) -7 (7+1 k 2 ) (r 2 +2fc 2 2 
(7-i fci) (7-1 fe 2 ) (iT+v^fci-v^fe) (ir+v^fei+v^fa) 

■y+i k 2 
7— i fc 2 

-2 7fe 1 2 +2ifci 3 +ifci (r 2 -2fc 2 2 )+ 7 (r 2 +2Ai2 2 ) 
(7-i fci) (ir+v^fei-v^fe) (iT+v^fci +v / 2fc 2 ) 

r 

r-i v / 2(fci-fc 2 ) 

r (T+i fci) (7+i fc 2 ) 

(7-1 fci) (7-1 fc 2 ) (r-iv^fci-iv 7 ^^) 

r(7+tfci) 

(7-jfci) (r-iv^fci+iv^fe) 
r( 7 +i k 2 ) 

(l-ik 2 ) (r-i^(ki+k 2 )) 

iV2(k 1 -k 2 ) 
T-iV2(k 1 -k 2 ) 

v / 2( 7 +»fc 2 ) (fcl+fcg) 
(7-« fa) (iT+v^fci+v^fe) 
V2(7+t fc 2 ) [fa+fcgj 
(7-1 fc 2 ) (iT+y^fci+v^fa) 
i y/2(fci-fc 2 ) 



(26) 



r-i-v/2(fci-fca) 

2 y^rfci (-i 7 +fc 2 ) 

(7-i fci) (iT+V2k!-V2k 2 ) (iT+V2 ki+V2 fc 2 ) 

2iy / 2"rfci (7+ifc 2 ) 

(7-i fci) (iT+v^fci-v^fa) (iT+v^fci+v^fa) 





At At = 0; B 5 = 1 we have: 



d r 

1 r-iV2k 1 +iy / 2k 2 

B = ~ iF 

2 iF+V2k ± +V2k 2 
r> _ it 

3 jr+ v / 2fci+v / 2fc 2 



B, 



6 

B 7 

B 8 
A 2 
A 3 
A, 



T-iVZk 






2 7 fci 2 + 


■2ifci 3 +ifci (r 2 -2fe 2 2 )-7(r 2 4 


-2fc 2 2 ) 


(7— i ki) 1 


^r+V2fci-V2fc 2 ) (ir+^2A,-iH 


-V2k 2 ) 


-2 7 fci 2 - 


-2ifci 3 -ifci (r 2 -2fc 2 2 )+7 (r 2 +2fc 2 2 ) 


(7— ifci) 1 


(irWIfci— \/2fc 2 ) (iT+^fciH 


-V2k 2 ) 


= -1 








2 v^Tfci (-i7+fc 2 ) 




(7— ifci) 1 


f AT+V2k 1 -V2k 2 ) (iT+v^fcH 


-V2k 2 ) 


= 








2 v^Tfci (i7+fc 2 ) 





(27) 



(7-ifci) (iT+v^fci-v^fe) (jT+v^fci+v^fe) 



_ ^(7+^ fcl)("fcl+fc 2 ) 



iT+v(2fci— v/2fc 2 

6 ~~ (7-jfci) (ir+ v /2fci+v / 2fc 2 ) 
J _ y/2(7+tfci) (fci-fc 2 ) 

7 (7-ifci) (ir+%/2fei— s/2fc2) 



V / 2(fcl+fc2) 



ir+v^fci+v 7 ^^ 



We are not able now to discuss the physical meaning of the results obtained 
in (j2H| and (J27)l because we need the similar expressions for the v + as part of the 
scattering amplitude. This is what we plan to do separately. However, some basic 
features of solutions obtained with (j2Tlj) and (|27j) can be already discussed. Here 
we shall pay attention to the scattering of the particles being initially bound in a 
cluster. The incoming cluster is represented by the incoming surface wave: 

e -ikx- ly . X = k 2_^2^ 

Substituting this into (|2l)j) one has at At = 1 and B 5 = 



Bt = l 



r 



i v / 2fc+v / 2 7-r 

B 2 = B 3 = 



B 4 = -1 



r 



-j V2(fc-« 7 )+r 



B 6 = 5 7 = 5 8 = 

. rt . _ . (k-ij) ('2fc 2 +4ifc7-27 2 +r 2 ) 

A 2 = A 3 = A 4 = 4 _ f; V ^-=1 — ^ ^ f98) 

2 d 4 (fc+i 7 ) (V2fc-iv / 27+ir)(v/2fc+j(v / 27+r)) l ZC V 

4r = F 

5 i v / 2fc+v / 27-r 

A 6 = 

^4 (»fc+ 7 )r 

7 (-jfc+7) (-i v / 2fc-v / 2 _ 7+r) 

A s = 



Considering (f2*U|) we should first mention that the system in question does not 
allow rearrangements, i.e. as the coupled particles come to the barrier they scatter 



back also coupled. No transitions like: 2 «-> 1+1; 2 <-> 1+1* is possible. Here * stands 
for the particle bound at the barrier. This exceptional, in general, but common for 
the Bethe ansatz situation is due to the symmetry restrictions we applied. The 
coefficient A4 which plays the role of the reflection coefficient depends effectively 
on the strength of both interaction with the barrier and the interaction between 
particles. However in order to study the influence of both interactions on the cluster 
conductance through the barrier one should combine the data obtained for t>_ with 
ones for v + . In the next section we shall start the solution of the scattering problem 
for v + . 

5 Symmetric case and the Riemann-Hilbert Prob- 
lem 

The symmetric case, i.e. the case which is center-mass symmetric with respect to 
the barrier, leads to the following expression 



through the 2x2 matrixes Ni(z) and N 2 (z) with polynomial coefficients and poly- 
nomials Pi(z) and P2(z). We do not present here the exact expressions for Ni(z) 
and Pi(z) we obtained by using the computerized tool for analytical calculations 
(Mathematica) . Instead of this we shall adopt our general results [HUH] to the case 
of the functional equation with 2x2 rational matrix coefficient as it has been done 
in (J29|) . In this way we can not obtain the exact solution for v + as it was done for 
f_, however we are able to study the properties of this solution in much details. 
Namely, for the application to the solution of scattering problem one needs to find 
analytic solutions of functional equations (fTSj) with given singularities, which would 
correspond to the incoming plane waves fUJ In order to do this it is useful to 
find a generating solution, i.e. a meromorphic matrix function having meromorphic 
inverse one with finite number of singularities in the band 9fcc G [0,1). Then the 
solution with given properties can be obtained through the action of generating 
solution on some periodic vector function which is simply a rational function of ex- 
ponentials. A generating solution for the functional equations (JT3j) can be obtained 
from the solution of the iterated equations (JUjl using the results of Theorem 1. So 
in the rest of the paper we shall concentrate on the functional equations (J 14)) only. 

We reformulate the equations (|14j) to the form of the boundary value problem 
by introducing the variable z = e 2 ^ , which coincide with the definition (|22|) in the 
case period 2n. 

Due to (|15|) the matrix function N (z) will be meromorphic, which is just obvious 
in our case, and the boundary values of solutions for (fT4"j) on the cut z G [0, +00) 
will satisfy 



This is the special case of the problem known as the matrix Riemann-Hilbert problem 




(29) 



g(z-iO) = N (z) g(z + iO) + E (z) . 



(30) 



. We propose the way to solve this type of boundary value problems and demonstrate 
the functional class of solutions. 

Let us define C := [0, +00), zi- := 0, zi + := +00, hi£ z := In z. Let us also note 
that the requirement N (zi±) = I is obviously fulfilled in our case. Here we need the 
results obtained for the Homogeneous Riemann-Hilbert Problem: Find the 
vector function g (z) |C 1— > C" analytic on C\ {£}, extendable until the both sides of 
C and such that its boundary values satisfy 

g(z)\ c+ = N(z)g(z)\ c _. (31) 

The solution of an inhomogeneous problem would be needed if we study some sin- 
gular solutions for v + (compare with P). 

In our consideration we use the matrix coefficient N (z) in the diagonal (or 
Jordan) form. Here we can omit the n x n case, which we briefly considered in our 
previous paper P and come directly to the simple diagonalization of a 2 x 2 matrix 
N. Below we shall use the symbol N(z) to refer any of A^/pi or N 2 /p 2 - Since any 
2x2 matrix has just two eigenvalues, which are given by the values of the two-sheet 
analytic function 

1 



A (z) = - [TrN (z) + ^JTy 2 N (z) - 4DetiV [z] 

the procedure of finding the diagonal or Jordan form is direct. Moreover, for the 
problem in question we can find explicitly the expressions for 



and 



TTiVi(z) — 4DetiVi(z) = 

-i28^ 4 (-i + 2 4 ) 2 7 2 r 2 

(1 + z 16 + (z 4 + z 12 )(- 7 4 + r 2 (4 + r 2 ) - 4 7 2 (i + 2r 2 ))- 
^ 8 (4 7 2 (2 + r 4 ) + 2(1 + 4r 2 + r 4 ) + 7 4 (2 + 4r 2 + r 4 ))) 



Tt 2 N 2 (z) - ADetN 2 (z) 
128z 4 (1 + z 4 ) 2 7 2 r 2 



(32) 



(33) 



(1 + z 16 + (z 4 + z 12 )( 7 4 - r 2 (4 + r 2 ) + 4 7 2 (i + 2r 2 ))+ 
z 8 (4 7 2 (2 + r 4 ) + 2(1 + 4r 2 + r 4 ) + 7 4 (2 + 4r 2 + r 4 ))). 

We conclude that we deal now with the case (l.b) of PP, Let us mention that the 
eigenvalues Ai (z) t^A 2 (z) a.e. and have branching points. 

The corresponding eigenprojections Aj (z) can be expressed as A, (z) = v» (z) x 
v* (z) where Vj (z) is the eigenvector of N (z) corresponding to the eigenvalue Aj (z) 
and \\ (z) is one for (z) . These eigenvectors are normalized as v* (z) x Vj (z) = 1, 
so TrAj (z) = 1, A 2 (z) = A { (z) and A^ (z) A { (z) = for % ^ j. 

Let us recall the notion of the generating solution we used pQ . We call an analytic 
matrix function X (z) generating solution of the Riemann-Hilbert problem 1 (|31|) 
if 



1 Note that the generating solution for the functional equation (|14|) can be obtained as X (e 27 ™ « ) . 



1. it satisfies the boundary value problem 

X(z)\ c+ = N(z)X(z)\ c _. (34) 

2. X (z) and X~ x (z) are regular on C/{£} except may be some finite number s 
of points {zi} s i=1 where they poles of finite multiplicities. 

Suppose we know some generating solution of the problem (|3ip. One can show 
that for any regular solution g (z) the function g (z) = X^ 1 (z) g (z) should satisfy 

go (z)\ c+ = go {z)\ c _ 

and so would be rational function with poles at the points {zi} only and of multiplic- 
ity not higher then one of X -1 (z). In general, any such function can be expanded 
in finite linear combination of simplest fractions and the regularity conditions at the 
points {z^ would lead to the homogeneous system of linear equations for the coef- 
ficients in this decomposition. However, we shall demonstration this fact elsewhere. 

In we proposed the way how to construct a generating solution for the general 
case of meromorphic matrix N(z). Here we shall show how this can be applied to 
the specific problem in question. To find the generating solution we make use of the 
fact that the matrix coefficient can be represented as 

N (z) = A+ (z) A + (z) + A_ (z) A_ (z) 

where the signs ± denote the values of functions on different sheets of the Riemann 
surface. In the considered case the structure of the Riemann surface is determined 
in or (jSni)? i- e - by the algebraic curves 



K x (z) = (1 + z ie + (z* + z 12 )(- 7 4 + r 2 (4 + r 2 ) - 4 7 2 (1 + 2r 2 ))- 

2 8 (4 7 2 (2 + r 4 ) + 2(1 + 4r 2 + r 4 ) + 7 4 (2 + 4r 2 + r 4 )))^ 



(35) 



or 



(36) 



K 2 (z) = (1 + z 16 + {z A + ^ 12 )( 7 4 - T 2 (4 + r 2 ) + 4 7 2 (1 + 2r 2 )) + 

z 8 (4 7 2 (2 + r 4 ) + 2(1 + 4r 2 + r 4 ) + 7 4 (2 + 4r 2 + r 4 )))^ 

which have by 16 branching points {hi}^ =1 in all non-degenerated cases, i.e. r ^ 0, 
7 7^ 0. The points hi are given by simple but bulk expressions we shall leave outside 
this text. Let us split the set {hi}]^ in two parts {hi}^ =1 = {ui}® =1 U {vi} i=1 and 
fix the cuts S{ which run between the points Ui and Vi so that they do not cross 
the half-axis [0, oo). Because the construction below is the same for both functions 
Ki(z) and K 2 (z) we shall omit the indices and refer these functions as K(z). We 
denote by K± (z) two branches of the function K(z). The branches can be fixed by 
the condition 

lim K + (z)/z 8 = 1. 

The corresponding branches of A (z) and A (z) are denoted as A± (z) and A± (z). If 
we denote the points of the Riemann surface 1Z of K (z) by z = (z,j); z G C; j = ± 
these branches are defined as X± (z) = A ((z, ±)) ; A± (z) = A ((z, ±)). 



Let us define the function g (z) on TZ 

j K(z) r+- \n\ + (z>) ,\ 
9o{z)=eXP \^nJo K + (z<){z<-z) dZ ) (37) 

We have to mention that (|37|) gives a new representation of go( z ) i n comparing with 
go(z) from pQ. This form is more compact and seems to be more suitable for the 
actual calculations since it does not require the calculation of A (z) singularities. 

The function go (z) is meromorphic on TZ\ {00} and has the essential singularity 
at z — ► 00. As z — > 00 the function g ( z ) behaves as e 9 ^ where the coefficients of 
the polynomial q (z) of the degree 8. 

This singularity can be canceled if one multiplies go ( z ) by a product of the Baker- 
Ahiezer (BA) function ip (z) of the Riemann surface TZ corresponding to the point 
00, the polynomial q(z) and some non-special divisor D of degree n Really, if 
one fixes n points on TZ and allows the function ip (z) to have simple poles at this 
points, i.e. specifies the divisor D, then the condition ip (z) ~ e~ q ^ at z ^00 will 
specify the function ip uniquely up to the constant factor. In order to get the explicit 
expression for the BA function one chooses the canonical sections {aj}^ =1 and {h}^ =1 
on the Riemann surface TZ. Then one constructs the second kind differential Q with 
the main part — dq (z) at 00 and normalized by the conditions 



n = o. 

When one denotes by U = {Ui, .., U m } the vector of b -periods of Q, 

Jbi 

and chooses some point zo on TZ the explicit expression for ip looks like ^Hl 

f r „\ (A(z) - A(D) + U - K) 

where K is the vector of Riemann constants, A denotes the Abel transform and 9 
is the Riemann theta-function of the surface TZ. 

Let us introduce the function g (z) = g Q (z) if) (z). The branches of g (z) on the 
sheets TZ we shall denote by 

9+ (z) = g((z,+)) ! <7_(z) = </((*,-)). 

In this notations the generating solution of the Riemann-Hilbert problem can be 
found as 

X(z)=g + (z) A + (z)+g- (z) A_ (z) . 

The regularization of go (z) described above has some shortcomings. Namely, the 
investigation of X~ l (z) meets the problem of theta-function zeroes and the formulae 
obtained are not a good base for the numerical calculations. In the concluding 
section of PQ we proposed another method to construct the BA function. The 
method leads to the set of Dubrovin's type [TH] ordinary differential equations. We 
plan, in the further study on the subject, to investigate the solutions of this equations 
for the particular case we study here. However we should emphasize that the method 
of the solutions construction for (fTU|) and their functional class is already studied. 



6 Conclusion 



In summary, we have studied a new mathematical model for correlated tunneling 
of electrons through a barrier in ID quantum wires. For one subproblem the exact 
solution has been obtained and for another one the reduction to the generalized 
Dubrovin's equation has been constructed. The latter can be solved either numer- 
ically or analytically. The future work will be devoted to the exact integration of 
the Dubrovin's equations obtained and to the calculation in terms of their solutions 
of the physical characteristics of correlated tunneling in quantum wires. 

Finally we note that the approach developed in this paper can be easily applied 
to the phenomena of type 1) mentioned in Introduction, i.e. to the study of exotic 
exciton complexes in ID semiconducting nanostructures. 
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